Solution of the center-of-mass problem in nuclear structure calculations 
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The coupled-cluster wave function factorizes to a very good approximation into a product of an 
intrinsic wave function and a Gaussian for the center-of-mass coordinate. The width of the Gaussian 
is in general not identical to the oscillator length of the underlying single-particle basis. The quality 
of the separation can be verified by a simple procedure. 
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The atomic nucleus is a self-bound quantum many- 
body problem. Its Hamiltonian is invariant under trans- 
lations and rotations. Thus, total momentum and to- 
tal angular momentum are conserved quantities, respec- 
tively. When one chooses a single-particle basis to solve 
the nuclear many-body problem, one is confronted with a 
dilemma. There is no single-particle basis of states that 
are simultaneous eigenstates of the momentum operator 
and the angular momentum operator. For calculations 
of finite nuclei within the nuclear shell model, one usu- 
ally chooses a spherical single-particle basis. Such a basis 
obviously breaks translational invariance, and the conse- 
quences have been addressed in numerous papers (see, 

e.g. na a is, 

Only two rigorous solutions to the problem are known. 
The first approach deals with the problem in translation- 
ally invariant Jacobi coordinates. We refer the reader to 
Refs @, a S for recent applications. Due to the factorial 
scaling of the required antisymmetrization, this method 
is restricted to few-body problems with A < 7 or so. 
The second approach is suitable for many-body problems 
but severely limits the choice of the single-particle basis 
to the eigenstates of the harmonic oscillator. For the 
A-body system, one considers an A-fermion basis that 
consists of all single-particle product states with exci- 
tations up to and including Nftw. The corresponding 
model space is called a full Nfuv space. In this space, 
eigenstates of the translationally invariant nuclear Hamil- 
tonian are also eigenstates of the center-of-mass Hamil- 
tonian H cm (ui) (see Eq.([3} below for a definition). Note 
that the Fock-space basis employed in this approach is 
not translationally invariant, i.e. the total momentum 
is not a conserved quantity. Rather, the complete Nhuj 
space ensures that all eigenstates are products of an in- 
trinsic state f/iin and a center-of-mass state ip C m- The in- 
trinsic states are invariant under translations, while the 
wave function of the center-of-mass coordinate is obvi- 
ously not an eigenstate of the total momentum. Thus, 
the factorization 



1p = V'cm^ii 



(1) 



is central to the Nhw space, and it is essential for a cor- 
rect treatment of the center-of-mass problem. Note that 



the particular form of the center-of-mass wave function 
tp cm is irrelevant [3,0] . Recent examples for this a ppr oach 
are, e.g., the no-core shell- model calculations rJ 11 1 and 
the Ohw shell-model calculations 13, l3 |. 

Unfortunately, the two rigorous approaches to the 
center-of-mass problem scale exponentially with the 
number of active nucleons. Furthermore, many nuclei of 
interest are very weakly bound, and the oscillator basis 
does not provide the correct radial asymptotics. Alterna- 
tive wave- function-based approaches such as the coupled- 
cluster method Q [H G3, |H GH El S, E| , the unitary- 
model operator approach [221 ] . or stochastic shell- model 
approaches 23|, 24 1 scale more gently with increasing size 
of the model space and the mass number. While these 
methods are also based on a single-particle basis, they 
do not employ a complete Nhu> space. Thus, there is no 
analytical guarantee that the wave functions computed 
by these methods also exhibit the factorization ((T|). For 
this reason, the results obtained by such methods are oc- 
casionally viewed with skepticism. The concern seems 
not primarily with the possible breaking of translational 
invariance, but rather with the problem to quantify the 
errors that might be involved. Recall that nuclear lattice 
calculations also break translational invariance, albeit in 
a controlled way pit ]. 

It is the purpose of this Letter to eliminate these con- 
cerns. Moreover, we present a simple tool (i.e., the com- 
putation of expectation values of the generalized center- 
of-mass Hamiltonian) that assesses to what degree the 
factorization ([T]) is exhibited in coupled-cluster calcula- 
tions. In what follows, we consider the 16 nucleus and 
demonstrate that the coupled-cluster wave function is a 
product of an intrinsic wave function and a Gaussian for 
the center-of-mass coordinate. 

For the nucleus 16 O we employ the low-momentum in- 
teraction Viowfe [26j], generated from Entem and Mach- 
leidt's chiral nucleon-nucleon potential [27] ] by imposing a 
smooth 28j momentum cutoff A = 1.8 fm _1 . Our single- 
particle basis consists of spherical oscillator states. The 
oscillator spacing fku and the size of the model space (in 
terms of the number N of major oscillator shells) are the 
parameters of our calculations. We employ the intrinsic 
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Hamiltonian 



T - T r , 



E 



l<i<j<A 



V . 

(pi 



2mA 



+ V(n-fj) 



(2) 



Here, T and V denote the kinetic and potential en- 
ergy operators, respectively, and T cm denotes the kinetic 
energy of the center of mass. The intrinsic Hamilto- 
nian is manifestly invariant under translations. We per- 
form a spherical Hartree-Fock calculation to obtain an 
optimized single-particle basis. This is followed by a 
spherical coupled-cluster calculation 29]. Clearly, the 
Hartree-Fock single-particle basis in combination with 
the coupled-cluster method in its singles and doubles ap- 
proximation (CCSD) does not employ a full Nhu space. 
Thus, the separation ([I]) is not guaranteed from the out- 
set of our calculations. 

Figure [T] shows the ground-state energy as a function of 
the oscillator spacing for a model space consisting of N = 
9 oscillator shells. The energy varies by about 1 MeV 
while the oscillator spacing hu varies by more than a 
factor of two. This shows that the ground-state energy 
is very well converged with respect to the size of the 
model space. Note that we could easily employ larger 
model spaces and improve the convergence of our results. 
However, this is not necessary for the purpose of this 
Letter, and we refer the reader to established benchmark 
calculations 130(1 . 
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FIG. 1: (Color online) Ground-state energy (within CCSD) 
of 16 O with a low- momentum potential as a function of the 
oscillator spacing hu). The model space consists of nine ma- 
jor oscillator shells. Inset: Expectation value E cln (u)) of the 
center-of mass Hamiltonian with the standard frequency de- 
pendence. 

Let us consider the generalized center-of-mass Hamil- 
tonian 



H cm (ui) — T a 



1 



mAu 2 R 2 n „ 



-hu 



(3) 



Here, Co is a free parameter and not necessarily identi- 
cal to the frequency u> of the underlying oscillator basis. 



The generalized center-of-mass Hamiltonian (|3]) exhibits 
a zero-energy Gaussian ground-state wave function for 
all values of a). In what follows, we demonstrate that the 
coupled-cluster ground state is the zero-energy eigenstate 
of H cm (Co) for a suitably chosen frequency Q. Thus, the 
coupled-cluster wave function factorizes, and the center- 
of-mass wave function is a Gaussian. 

We denote the expectation value of the generalized 
center-of-mass Hamiltonian ([3]) in the coupled-cluster 
ground state as 



E cm (uj) = (H cm (ui)) 



(4) 



We compute the expectation value ((4| via the Hellmann- 
Feynman theorem, i.e. we add a small perturbation 
(3H cm (u) with — 0.001 to the intrinsic Hamiltonian ((2]), 
and the expectation value results from the difference quo- 
tient U, 32j ]. 

First, we consider the standard center-of-mass Hamil- 
tonian and set to = u>. The results are shown in the 
inset of Fig. [TJ Comparing the results for the center-of- 
mass energy E cm (u) with the ground-state energy sug- 
gests that the latter does not depend on the former. This 
is a first hint that the intrinsic and center-of-mass coordi- 
nates decouple. For the model space with hu m 20 MeV 
the expectation value E cm (u) vanishes approximately, 
thus indicating that the wave function factorizes. As 
a check, we fix hu = 20 MeV, add the term f3H cm (u) 
to the intrinsic Hamiltonian, and compute the ground- 
state energy of the resulting Hamiltonian. We find that 
the ground-state energy varies by merely 15 keV as f3 
is increased from zero to one. Since the intrinsic en- 
ergy shown in Fig. [T] is practically independent of the 
frequency of the underlying oscillator basis, the wave 
function must approximately factorize for all values of 
hu shown in Fig. [T] However, only for hu w 20 MeV, 
do we know that the center-of-mass wave function is a 
Gaussian. 

Let us assume that the center-of-mass wave function 
generally has a Gaussian shape, i.e. it is the ground 
state of the generalized center-of-mass Hamiltonian Q 
for a suitably chosen frequency Co. It is thus our task 
to determine this frequency. To this purpose, we employ 
the identity 



H cnl (u) + \hu - T cm = ^ 
2 u l 



i(^) + -fi^ - T C] 



take its expectation value, require E cm (u) = 0, employ 
(T cm ) = § Hu), and determine the unknown frequency u 
from the already computed expectation values E cm (u). 
This yields the two possible frequencies 



2 /4 
hu = hu + -E cm (u) ± W -(E c 



t (w)) 2 + -huE cm (u) .(5) 



We employ these frequencies in the generalized center- 
of-mass Hamiltonian and compute the corresponding ex- 
pectation values E cra (u). We find that one expectation 
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value is close to zero, while the other is usually very large. 
The small expectation values are shown in Fig. [2 and 
the corresponding frequencies Co are shown in the inset of 
Fig. O The expectation values are very small compared 
to the energy tub « 20 MeV of spurious center-of-mass 
excitations. The practically vanishing expectation values 
demonstrate that the coupled-cluster wave function fac- 
torizes. This is the main result of this Letter. The reader 
should not be concerned about the fact that some of the 
expectation values shown in Fig. [2] assume small negative 
values (of size -0.01 MeV). Recall that the coupled-cluster 
method is non-variational when the cluster operator is 
truncated, and that small negative expectation values are 
thus tolerable. As shown in the inset of Fig. [2j the fre- 
quency Co corresponding to the Gaussian center-of-mass 
wave function varies only little as the frequency u> of the 
underlying oscillator basis is changed. We made the fol- 
lowing two checks. First, we computed the expectation 
value of T cm and found that T cm « fto, as expected for a 
Gaussian. Second, we repeated the calculations directly 
in the oscillator basis and did not employ the Hartree- 
Fock basis. We again find very small values for E cm (Co) 
of order -0.1 MeV, and an almost constant frequency Co 
very close to what we found in the Hartree-Fock basis. 
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FIG. 2: (Color online) 1 O ground-state expectation value 
Bcm(ii) (within CCSD) of the generalized center-of-mass 
Hamiltonian H C m(u>) as a function of the oscillator spacing 
hui. The model space consists of nine major oscillator shells. 
Inset: Relation between the frequency ui and the frequency uj 
of the underlying oscillator basis. 



Let us also consider interactions with higher momen- 
tum cutoff. Wc employ Entem and Machleidt's chiral 
nucleon-nucleon interaction 27| at next-to-next-to-next- 
to-leading order (N 3 LO). This interaction has an approx- 
imate high-momentum cutoff of A m 500 MeV. We apply 
the coupled-cluster method to compute the ground-state 
energy of 16 O. Due to the relatively high momentum cut- 
off of the interaction, the wave function is more corre- 
lated. This requires us to employ a large model space 
and three-particle-three-hole cluster amplitudes to ob- 
tain converged solutions. The model spaces consist of 



N = 19 oscillator shells; however, the maximal single- 
particle orbital angular momentum is kept at I < 13. We 
have verified that this is sufficient for a convergence of 
the results. For the three-particle-three-hole clusters, we 
employ the A-CCSD(T) approximation [33, 34 1 . 

The bottom part of Fig. [3] shows that the ground-state 
energy is practically independent of the oscillator fre- 
quency of the underlying single-particle basis for a large 
frequency range. This demonstrates that the results are 
well converged with respect to the size of the model space. 
We find that the expectation value E cm (ui) of the stan- 
dard center-of-mass Hamiltonian increases with increas- 
ing frequency to of the model space, and that it assumes 
values of tens of MeV and varies strongly with the fre- 
quency to of the underlying oscillator basis. Clearly, the 
intrinsic ground-state energy is independent of the expec- 
tation value E cm (u>), and this suggests again a decoupling 
of the intrinsic and the center-of-mass wave functions. 

Let us again assume a Gaussian shape for the center- 
of-mass wave function and follow our two-step procedure. 
First, we employ Eq. §5§ and compute the two possible 
frequencies Co that are consistent with the already com- 
puted expectation value E cm (u>). Second, we compute 
the two corresponding expectation values E cla (Co). As 
before, we find a large and a small expectation value for 
each frequency to of the underlying oscillator basis. The 
small expectation values are shown in the top part of 
Fig. [3] The corresponding frequencies Co are shown in 
the middle part of Fig. [3l While not as impressive as for 
the low-momentum interaction, the expectation values 
E C m(Co) are below 1 MeV in size, and much smaller than 
the binding energy of 16 O or any of its excitations. In 
particular, the expectation values are small compared to 
the energy hCo rj 16 MeV of the spurious center-of-mass 
excitations. A simple two-level model yields that the 
wave function has about 6% squared overlap with spu- 
rious states. The relative small negative eigenvalues are 
again due to the non-variational character of the coupled- 
cluster method, and they are tolerable within the overall 
accuracy of the calculation. Thus, we can conclude that 
the coupled-cluster wave function also exhibits an ap- 
proximate factorization ([1]) for the interaction with a rela- 
tively high momentum cutoff. We speculate that quadru- 
ple cluster amplitudes would be necessary to further re- 
duce the magnitude of E cm (Co). Note that the frequency 
Co that determines the width of the Gaussian center-of- 
mass wave function is again almost constant over a wide 
range of frequencies of the underlying oscillator basis. 

Finally, we employed the chiral potential also for 
4 He and show the results in Fig. 0J The bottom part 
shows the ground-state energy (E « —25.56 MeV in A- 
CCSD(T)) and compares it to virtually exact Faddeev- 
Yakubowski calculations (E « -25.41 MeV [35]). The 
middle part shows the frequency of the approximate 
Gaussian center-of-mass wave function, while the top 
part shows the expectation value E cm (Co). Again, the 
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FIG. 3: (Color online) Bottom: ground-state energy of 16 O 
within the A-CCSD(T) approximation as a function of the 
frequency hui of the underlying oscillator basis. Middle: Re- 
lation between the frequency ui and the frequency u> of the 
underlying oscillator basis. Top: Expectation value E cul (ui) 
of the center-of-mass vs. fou). 



approximate factorization is very statisfactory. At this 
moment, we have no profound understanding of the ob- 
served factorization. 
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FIG. 4: (Color online) Same as Fig. E] but for 4 He. 

In summary, we presented strong numerical evidence 
that wave function of the atomic nucleus is a product of 
an intrinsic and a center-of-mass wave function in a suf- 
ficiently large model space even when the model space 
is not a complete Nhui oscillator space. The center-of- 
mass wave function is approximately a Gaussian whose 
width varies little with the frequency u) of the underlying 
oscillator basis. The reported results open the door for 
a verifiable description of translationally invariant states 
for a large variety of model spaces and many-body meth- 
ods. 
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